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Abstract 

The effective perturbation theory developed by Braaten and Pisarski for 
gauge theories at finite temperature is extended to finite chemical potential. 
As a first application the collisional energy loss of a heavy quark propagating 
through a quark-gluon plasma with non-vanishing quark chemical potential is 
considered. Assuming /_i/T ~ 1, motivated by numerical simulations of heavy 
ion collisions at RHIC energies, we find that the effect of the quark chemical 
potential is rather small, unless the energy density instead of the temperature 
is fixed. 
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I. INTRODUCTION 



QCD lattice calculations predict a phase transition from hadronic matter to a phase 
of deconfined quarks and gluons (QGP) above a critical energy density Cc ~ 1 GeV/fm^ 
corresponding to a critical temperature Tc ~ 150 MeV (see, e.g. The creation of a QGP 
is expected to happen in ultrarelativistic heavy ion collision and in the early universe shortly 
after the big bang. An overview of the present status of the experiments is given in |0. 

Depending on the available center of mass energy of the colliding nuclei, one can think 
of two scenarios |Q: Below an energy per nucleon of approximately 5 GeV/A complete 
stopping of the participating nuclei takes place. It is believed that in this region of high 
compression and high temperature a QGP containing more quarks than anti-quarks might be 
produced. Above an energy of about a few 100 GeV/A, on the other hand, full transparency 
is expected. The QGP is believed to be created in the region between the two expelling and 
highly excited nuclei by vacuum polarization leading to a vanishing baryon number in the 
fireball. 

We consider the quark chemical potential /i, which is related to the baryo-chemical 
potential fit through fi = ^b/^ and measures the deviation from the balance of quarks and 
anti-quarks. An equal amount of quarks and anti-quarks corresponds to /i = 0, whereas 
^ > describes an excess of quarks over anti-quarks. It is generally agreed that /i/T ~ at 
LHC energies. At RHIC energies, however, there might still exist a considerable amount of 
stopping leading to fi/T ~ 1 - 2 0] as predicted by numerical simulations (RQMD) p. 

Properties of a QGP at finite chemical potential cannot be described by lattice calcu- 
lations so far. Perturbative QCD at finite temperature and chemical potential suffers from 
infrared singularities and gauge dependent results. These problems can be avoided for the 
most physical quantities of interest by adopting an effective perturbation theory developed 
by Braaten and Pisarski, which is based on the use of resumed Green functions instead of 
bare ones 0. In this way consistent results, i.e. complete to leading order in the coupling 
constant and gauge independent, can be derived in the weak coupling limit, g <^ 1. 
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Here we will consider the Braaten-Pisarski method at finite temperature as well as chem- 
ical potential, which is relevant for ultrarelativistic heavy ion collisions if no complete trans- 
parency occurs. Up to now only the cases T 7^ 0, /i = and T = 0, /i 7^ in an 
astrophysical context |j^ have been discussed. 

The next section treats the technical aspects of introducing the chemical potential into 
thermal field theory. We extend the Braaten-Pisarski method to finite chemical potential 
and apply the results to the calculation of the energy loss of a heavy quark in a QGP as 
a typical example in the third section. Throughout the paper we use the metric tensor 

g = {-{ ), natural units {h = c = k = 1) and the notation K = [ko, k) and = |k| for 

the momenta. 

II. PERTURBATIVE TECHNIQUES AT FINITE CHEMICAL POTENTIAL 

A. Saclay method 

The Braaten-Pisarski method starts from the imaginary time formalism where one has to 
sum over discrete energies (Matsubara frequencies). For this purpose the Saclay method was 
developed in the sixties by Balian and De Dominies and used by Pisarski in allowing 
an elegant way to evaluate finite temperature loops. The idea is based upon the 'mixed' 
representation of the propagators, i.e. a momentum representation of space components and 
a space representation of the time component. In the case of the bare fermion propagator 
including the chemical potential and neglecting the bare mass it reads: 



S{K) 
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The sum is evaluated best by a complex integration [|T0| yielding the Saclay representation 
of the fermion propagator: 

A(r, k) = ^ {e-'^ [1 - riFik - /x)] - e'^ nj.{k + /x)} , (2) 

where np{k ± /x) = l/[exp((A; ± /u)/T) + 1] denotes the Fermi distribution function. The 
advantage of this representation is the easy summation over the Matsubara frequencies since 
the discrete energy k^ appears only in exp(A;oT) in the propagator A(i^) . The derivation of 
the boson propagator containing no chemical potential proceeds analogously 



B. Hard Thermal Loops 

Following the program given by the Braaten-Pisarski resummation technique in ref. 
we calculate the Hard Thermal Loops (HTL) at finite chemical potential, first. HTL are 
one-loop diagrams for which only the contribution from hard loop momenta of the order T 
or larger is considered. This approximation is equivalent to the leading term of the high 
temperature expansion of the diagram under consideration. Here, we attach importance 
mainly to changes introduced by the chemical potential and refer to the literature otherwise. 

First we discuss the HTL approximation of the gluon self energy, including the quark 
chemical potential and using the Saclay method. The only modification to the HTL po- 
larization tensor compared to the zero chemical potential case comes from the diagram of 
fig.l. Using standard Feynman rules in the imaginary time formalism and the Saclay repre- 
sentation of the propagators, where we consider one fermion line as a quark with chemical 
potential /i and the other one as an anti-quark with — /i, we obtain analogously to the case 
of zero chemical potential assuming the HTL limit. 



HL(ir)=HooW 



ko , fco + 



_3^2 (3) 

3 \ 2k ko-k ^ ^ 



for the longitudinal part and 
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o 2 '^0 



(4) 



for the transverse part of the gluon polarization tensor. Here 

'Vz-iVz + lEfl. (5) 

is called the effective gluon mass, where Nf denotes the number of thermal flavors in the 
QGP and fif the corresponding chemical potentials. In the following a QGP consisting of 
up and down quarks with yU« = /id = /i and all other /ij = is assumed. The only change in 
the gluon self energy due to the introduction of a chemical potential is given by the effective 
gluon mass which depends now on the chemical potential in accordance with the result of 
the high temperature limit [|lT|. In the limit T ^ the result (§) to (^) also reduces to the 



one found in ref. 0. 

The same behavior is observed for the HTL quark self energy. One simply has to replace 
the effective quark mass by the chemical potential dependent one 



2 g^T^ f 1 fi\ 



already quoted in ref. |T2|. (In the course of the derivation of the HTL quark self energy 
we had to assume that the chemical potential is not much larger than the temperature, 
/i ^ T/g, whereas this restriction was not necessary for the gluon self energy.) 

Owing to Ward identities the HTL self energies also imply the existence of HTL vertex 
corrections @]. They also will be modified in the presence of a finite chemical potential in a 
similar way as the self energies. However, for the most applications considered so far they 
can be neglected. 
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C. Effective propagators 



The effective gluon and quark propagators follow by resuming the HTL self energies by 
means of a Dyson- Schwinger equation B. Hence they differ from the ones at zero chemical 



potential (see, e.g. ref. |13| for the gluon propagator and ref. |14] for the quark propagator) 
only by the effective masses (|^) and (^. According to the Braaten-Pisarski method effective 
Green functions (propagators or vertices) have to be used only if all legs of the Green function 
have soft momenta of the order gT. Otherwise the use of bare Green functions is sufficient. 
In this way gauge independent results, which are complete to leading order in the coupling 
constant, are obtained for physical quantities. Moreover, due to the self energies in the 
denominator of the effective propagators important medium effects such as Debye screening, 
which improves the infrared behavior, and Landau damping are included. Thus the Braaten- 
Pisarski method means a crucial improvement compared to the naive perturbation theory 
based on bare Green functions. However, the application of the Braaten-Pisarski method is 



limited by two problems [0. First, the Braaten-Pisarski method relies on the distinction 
between hard and soft momenta, i.e. on the weak coupling limit assumption, g <^ 1. The 
validity of the extrapolation of the results to realistic values of the coupling constant, g = 1.5 
- 2.5, in ultrarelativistic heavy ion collisions, has to be checked for each quantity individually. 
In general it appears to work for particles with large momenta, p ^ T. Secondly, the effective 
gluon propagator suffers from the absence of a static magnetic screening mass in the HTL 
approximation of the gluon polarization tensor, thus still leading to infrared singularities in 
certain quantities, as e.g. the damping rate of a hard parton [Q. However, so far only a 



few observables, e.g. the color relaxation are known to be affected by this problem. 
The most physical quantities of interest, such as the collisional energy loss of energetic 
partons [|18[, transport rates of quarks and gluons [1^, the viscosity of the QGP [0, and 
the production rate of hard photons in the QGP |21[], are to leading order infrared finite by 



using effective propagators due to dynamical magnetic screening |22 



Summarizing this section, the generalization of the Braaten-Pisarski method to finite 
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chemical potential is straightforward by modifying the effective parton masses, containing 
now also the chemical potential, that enter into the effective propagators. 



III. ENERGY LOSS OF A HEAVY QUARK 

The energy loss of energetic partons in a QGP is related to the phenomenon of jet 
quenching in ultrarelativistic heavy ion collisions which might serve as a signature for 
the QGP formation [^]. The energy loss of a parton in a QGP is caused by elastic scattering 
off the thermal partons (coUisional energy loss) and gluon bremsstrahlung (radiative energy 
loss). The coUisional contribution follows from the damping rate by introducing the energy 
transfer per parton collision as a weight under the integral defining the damping rate . 



The damping rate of a hard quark follows from the imaginary part of the quark self 
energy diagram shown in fig. 2, which comes from the imaginary part of the effective gluon 
propagator corresponding to Landau damping of a virtual collective gluon mode. This damp- 
ing mechanism is equivalent to elastic scattering off the thermal partons via the exchange 
of a collective gluon (plasma wave), as can be seen from cutting the diagram in fig.2. As 
mentioned above, the damping rate, which is quadratic infrared divergent using bare prop- 
agators, turns out to be logarithmically infrared divergent applying the Braaten-Pisarski 
method reflecting the absence of static magnetic screening in the effective gluon propagator. 
Integrating over the loop momentum a cancellation of the effective gluon mass takes place 
||T9| . Consequently the damping rate of a hard quark, given to logarithmic accuracy by 
7g = (7^T/(37r) ln{l/g) p5[, does not depend on the chemical potential. (The Fermi-Dirac 
distributions from the internal quark line is exponentially suppressed for hard quarks [0.) 

The additional energy transfer factor in the loop integral reduces the infrared divergence 
and leads to a finite result using the effective gluon propagator at soft momentum transfer, 
q^gT. For large momentum transfer, g~T, the exchange of a bare gluon in the elastic 



scattering diagram is sufficient. According to the prescription by Braaten and Yuan |26 



the both contributions can be matched onto each other by assuming a separation scale, 
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gT ^ g* ^ T, in the weak coupling limit. Adding up the soft and hard contributions, 
the separation scale drops out of the final result. In this way a infrared finite and gauge 
independent result for the coUisional energy loss of a heavy quark is found, which is complete 
to leading order of the coupling constant and independent of arbitrary parameters such as 
infrared and ultraviolet momentum cut-offs or separation scales. 

Taking a finite quark chemical potential into account, the derivation of the soft part 
proceeds analogously to the case of vanishing chemical potential [|T3i|l8| . The only change 
comes from the effective gluon mass that does not cancel as opposed to the case of the 
damping rate. The hard part requires integrations over Fermi-Dirac distributions, which 
can be achieved only numerically at finite chemical potential. It involves integrals of the 
form 



where 7 = 0.57722... is Euler's constant, ({z) Riemann's zeta function with ('(2)/({2) = 
—0.56996..., and the function F{fi/T) with -F(O) = 0, shown in fig. 3, has to be determined 
numerically. 

The final result for the energy loss of a heavy quark with mass M ^ T, energy E, and 
velocity f in a QGP of Nf thermal flavors after adding up the soft and hard contributions 
is given by rather simple expressions in the limits E <^ M'^/T, 






and E^ M'^ 



IT. 



dE SnalT^ 



dx 3 
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respectively, where B{v) is a smooth function of the velocity that increases monotonically 
from -B(O) = 0.604 to a maximum of 5(0.88) = 0.731, and then decreases to B{1) = 0.629 
p!8| . The expressions (||) and (||) extend the Bethe-Bloch formula |2^ to the case of a 
relativistic plasma with finite temperature and chemical potential. They differ from the 
corresponding expressions at zero chemical potential |TB| by the effective gluon mass (^, 
the effective flavor number Nf, and the additional function F{fi/T). The details of the 
derivation of (H) and (||) are presented in ref. ^8 . 



Figs. 4 and 5 show the energy loss of a charm and bottom quark extrapolated to = 0.2 
at T = 250 MeV and Nf = 2 for fi/T = 0, 1, and 2 as a function of the momentum p. 
In the case of a charm quark is matched onto (|^) at some intermediate momentum by 
demanding dE/dx to be continuous, whereas for a bottom quark (P) is sufficient for the 
momentum range chosen in fig.5. At momenta below a few GeV the energy loss turns out to 
be negative reflecting the breakdown of the extrapolation to realistic values of the coupling 
constant for small momenta |T5| , p!8| . The enhancement of dE/dx with increasing chemical 
potential is rather small (below about 25% at /i/T = 2) since there is a partial cancellation 
between the factor Nf increasing with /i and the function F{fi/T) which is negative for 
/i > (fig. 3). This demonstrates that the temperature is more important than the chemical 
potential for the energy loss. 



However, if one fixes the energy density, given e.g. by the MIT-bag equation of state 29 



377r^ llvraA . „ / 2a;, \ o 3 2a 



with the bag constant B = (200 MeV^), increasing the chemical potential corresponds to a 
decrease of the temperature resulting in a significant reduction of the energy loss as shown 
in fig. 6. (At /i = a temperature of T = 250 MeV corresponds to e = 5.232 GeV/fm^.) 
Considering ultrarelativistic heavy ion collisions the choice of a certain energy density instead 
of a temperature might be reasonable as the former is closer related to the collision energy 



The collisional energy loss of a light quark follows from (^) by replacing the energy Ehj a 



maximum momentum transfer g ~ E/2, the energy loss of a gluon by multiplying the one of 
a light quark by a color factor of 9/4 [^]. The radiative energy loss by gluon bremsstrahlung 
is much more complicated and has not been addressed by the Braaten-Pisarski method so 
far. Estimates based on naive perturbation theory indicate the dominance of the radiative 
energy loss over the coUisional one in a QGP 



IV. CONCLUSIONS 

In the present work the Braaten-Pisarski method has been extended to finite chemical 
potential. The effective propagators, based on a resummation of the HTL self energies, differ 
from the one at zero chemical potential only by the effective masses which now depend on the 
chemical potential. The Braaten-Pisarski method can be used for calculating consistently 
quantities of the QGP that are logarithmically infrared divergent if bare Green Functions 
are utilized. Interesting observables related to the evolution and detection of a QGP in 
ultrarelativistic heavy ion collisions, belonging to this class of quantities, are the energy loss 
of energetic partons propagating through the QGP, the transport rates of quarks and gluons, 
the viscosity of the QGP, and the production rate of hard photons in the QGP. 

Here the coUisional energy loss of a heavy quark, representing a typical example for 
the application of the Braaten-Pisarski method, has been computed to leading order at 
finite temperature as well as finite quark chemical potential. This might be relevant for 
heavy ion collisions at RHIC energies, for which numerical simulations (RQMD) predict 
a considerable amount of stopping leading to a quark chemical potential of the order of 
/i ~ 1 - 2 T. The energy loss of a heavy quark has been calculated consistently adopting 
the Braaten-Yuan prescription, which corresponds to a decomposition in a soft contribution, 
calculable with an effective gluon propagator, and a hard, within naive perturbation theory 
accessible contribution. Besides the modification caused by the effective gluon mass 
in the soft part, momentum integrations over the Fermi-Dirac distribution functions in the 
hard part lead to an additional term (see (0)), which requires a numerical evaluation at finite 
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chemical potential. However, these changes are small, typical below 25%, for ^ < 2T as long 
as the temperature is fixed, showing the dominant role of the temperature over the chemical 
potential. If, on the other hand, a given energy density e{fi, T) is assumed, corresponding to 
a decrease of the temperature with increasing chemical potential, a significant suppression 
is observed. 

The same conclusions apply to the transport rates of quarks and gluons, which determine 
thermalization times and the viscosity of the QGP [|l9l. The transport rate follows from the 
damping rate by including a transport weight describing momentum relaxation. Owing to 
this transport weight the computation of the transport rate proceeds similarly to the one 
of the energy loss via the Braaten-Yuan prescription. The details of this calculation can be 
found in ref. p8| . 

The production rate of hard photons from a QGP can be derived to leading order again 
by using the Braaten-Yuan prescription, where the soft part follows from the photon self 
energy containing one effective quark propagator and the hard part from the tree matrix 
elements of Compton scattering and pair annihilation involving one gluon and one photon 
|2T|. The calculation of the hard photon production rate at finite quark chemical potential 



considering the effective quark propagator modified by the effective quark mass (|^) will be 
presented in a in a forthcoming publication |]33| . 
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FIGURES 

FIG. 1. Part of the HTL gluon polarization tensor containing the quark chemical potential 
FIG. 2. Quark self energy containing the effective gluon propagator 
FIG. 3. The function F(^/T) 
FIG. 4. Energy loss of a charm quark for /x/T = 0, 1 and 2 for a given temperature 

FIG. 5. Energy loss of a bottom quark for ji/T = 0, 1 and 2 for a given temperature 

FIG. 6. Energy loss of a bottom quark for ///T = 2 for a given energy density compared to a 
given temperature 
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